Limits on temporal variation of quark masses and strong interaction from atomic 

clock experiments 
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We perform calculations of the dependence of nuclear magnetic moments on quark masses and 
obtain limits on the variation of (mq/ A.qcd) from recent atomic clock experiments with hyperfine 
transitions in H, Rb, Cs, Elg^ and optical transtion in Elg^. 
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I. INTRODUCTION 

Interest in the temporal and spatial variation of major 
constants of physics has been recently revived by astro- 
nomical data which seem to suggest a variation of the 
electromagnetic constant a — /Tic at the 10~^ level for 
the time scale 10 billion years, see 0) (a discussion of 
other limits can be found in the review |3| and references 
therein). However, an independent experimental confir- 
mation is needed. 

The hypothetical unification of all interactions implies 
that variation of the electromagnetic interaction constant 
a should be accompanied by the variation of masses and 
the strong interaction constant. Specific predictions need 
a model. For example, the grand unification model dis- 
cussed in predicts that the quantum chromodynamic 
(QCD) scale Aqcd (defined as the position of the Lan- 
dau pole in the logarithm for the running strong coupling 
constant) is modified as follows 



5A.QCD _ g^fe 



(1) 



The variation of quark and electron masses in this model 
is given by 



~ 70 — 

m a 



(2) 



This gives an estimate for the variation of the dimension- 
less ratio 



S{m/AQCD) ^ ^^Sa 
(m/AQCD) a 



(3) 



The large coefficients in these expressions are generic for 
grand unification models, in which modifications come 
from high energy scales: they appear because the running 
strong coupling constant and Higgs constants (related to 
mass) run faster than a. This means that if these models 
are correct the variation of masses and strong interaction 
may be easier to detect than the variation of a. 

Unlike for the electroweak forces, for the strong inter- 
action there is generally no direct relation between the 
coupling constants and observable quantities. Since one 



can measure only variation of the dimensionless quan- 
tities, we want to extract from the measurements vari- 
ation of the dimensionless ratio mq/ Aqcd where mq is 
the quark mass (with the dependence on the normaliza- 
tion point removed). A number of limits on variation 
of mq/ Aqcd have been obtained recently from consid- 
eration of Big Bang Nucleosynthesis, quasar absorption 
spectra and Oklo natural nuclear reactor which was ac- 
tive about 1.8 billion years ago d H H (see also 
11 EIH Il3)- Below we consider the limits which 
follow from laboratory atomic clock comparison. Lab- 
oratory limits with a time base about a year are espe- 
cially sensitive to oscillatory variation of fundamental 
constants. A number of relevant measurements have been 
performed already and even larger number have been 
started or planned. The increase in precision is very fast. 

It has been pointed out by Karshenboim that mea- 
surements of ratio of hyperfine structure intervals in dif- 
ferent atoms are sensitive to variation of nuclear mag- 
netic moments. First rough estimates of nuclear mag- 
netic moments dependence on mq/ Aqcd and limits on 
time variation of this ratio have been obtained in our 
paper Using H, Cs and IIg+ measurements [T^ [T5j . 
we obtained the limit on variation of mq/ Aqcd about 
5 • 10"^'^ per year. Below we calculate the dependence of 
nuclear magnetic moments on mq/ Aqcd and obtain the 
limits from recent atomic clock experiments with hyper- 
fine transitions in H, Rb, Cs, Hg"*" and optical transition 
in IIg+. It is convenient to assume that the strong inter- 
action scale Aqcd does not vary, so we will speak about 
variation of masses. 

The hyperfine structure constant can be presented in 
the following form 



A — const X 



\[a'^Frei{Za)][fi 



M„ 



(4) 



The factor in the first bracket is an atomic unit of en- 
ergy. The second "electromagnetic" bracket determines 
the dependence on a. An approximate expression for the 
relativistic correction factor (Casimir factor) for s-wave 
electron is the following 



rel 



7(47' - 1) 



(5) 



where 7 ~ y/l — (Za)^, Z is the nuclear charge. Varia- 



2 



tion of a leads to the following variation of Frei 0| : 



K 



= A 

Frel OL 



72(472 - 1) 



(6) 



(7) 



More accurate numerical many-body calculations 16J of 
the dependence of the hyperfine structure on a have 
shown that the coefficient K is slightly larger than that 
given by this formula. For Cs [Z—hh) K— 0.83 (in- 
stead of 0.74), for Rb K=^M (instead of 0.29), for Hg+ 
K=2M (instead of 2.18). 

The last bracket in eq. I0J contains the dimcnsionless 
nuclear magnetic moment /i in nuclear magnetons ( the 
nuclear magnetic moment M = M^ITc)' electron mass 
rrie and proton mass Mp. We may also include a small 
correction due to the finite nuclear size. However, its 
contribution is insignificant. 

Recent experiments measured time dependence of the 
ratios of hyperfine structure intervals of i^^Hg+ and H 
[T^ . ^-^^Cs and ^''Rb and ratio of optical frequency 
in Hg+ and ^■^^Cs hyperfine frequency [ig- In the ratio 
of two hyperfine structure constants for different atoms 
time dependence may appear from the ratio of the fac- 
tors Frel (depending on a) and ratio of nuclear magnetic 
moments (depending on niq/ Kqcd)- Magnetic moments 
in a single-particle approximation (one unpaired nucleon) 
are: 



li^{gs + (2j - l)gi)/2 



for 3^1 + 1/2. 



M = 



2(J + 1) 



{-Qs + (2j + 3)50 



(8) 



(9) 



for j = / — 1/2. Here the orbital g- factors are gi — 1 
for valence proton and gi ^ for valence neutron. The 
present values of spin g- factors gs are gp = 5.586 for 
proton and (?„ = —3.826 for neutron. They depend on 
Tfiq/ Aqcd- The light quark masses are only about 1% 
of the nucleon mass (m^ = (r7i„ -I- md)/2 « 5 MeV). 
The nucleon magnetic moment remains finite in the chiral 
limit of m„ — nid — 0. Therefore, one may think that the 
corrections to gs due to the finite quark masses are very 
small. However, there is a mechanism which enhances 
quark mass contribution: 7r-meson loop corrections to 
the nucleon magnetic moments which are proportional 
to TT-meson mass m^r ~ ^Jrr^KqcD] m^ = 140 MeV is not 
so small. 

According to calculation in Ref. dependence of the 
nucleon g-factors on 7r-meson mass ttItt can be approxi- 
mated by the following equation 



g{m^) 



5(0) 



bm'i 



(10) 



where a= 1.37/GeV , b= 0.452/GeV2 for proton and a= 
1.85/GeV , b= 0.271/GeV2 for neutron. This leads to 
the following estimate: 



Sgp 

gp 

Sgn 

gn 



-0.174- 



-0.213 



-0.087^ 



-0.107^ 



(11) 



(12) 



Eqs. I|8I9I11I12|) give variation of nuclear magnetic mo- 
ments. For hydrogen nucleus (proton) 



^ ^ ^ ^ _0 087^ 
gp ' niq ' 



(13) 



For '^^^Hg we have valence neutron (no orbital contribu- 
tion), therefore the result is 



M gn 



Sm„ 
-0.107 ^ 



(14) 



For ^^Cs we have valence proton with j—7/2, l=A and 



^^0.22^=0.11^ 

jl JTItt niq 



(15) 



For °'Rb we have valence proton with ^=3/2, l—l and 



^ = _0.128^ = -0.064^ 

fl TOtt mq 



(16) 



Deviation of the single-particle values of nuclear mag- 
netic moments from the measured values is about 30 %. 
Therefore, we tried to refine the single-particle estimates. 
If we neglect spin-orbit interaction the total spin of nu- 
cleons is conserved. The magnetic moment of nucleus 
changes due to the spin-spin interaction because valence 
proton transfers a part of its spin < > to core neu- 
trons (transfer of spin from the valence proton to core 
protons does not change the magnetic moment). In this 
approximation gs — {1~ b)gp + bgn for valence proton (or 
gs — {1 — b)gn + bgp for valence neutron). We can use 
coefficient 6 as a fitting parameter to reproduce nuclear 
magnetic moments exactly. The sign of gp and gn are 
opposite, therefore a small mixing 6 ~ 0.1 is enough to 
eliminate the deviation of the theoretical value from the 
experimental one. Note also that it follows from eqs. (jlll 



This produces an additional suppres- 



[Tgll that ^ 

sion of the effect of the mixing. This indicates that the 
actual accuracy of the single-particle approximation for 
the effect of the spin g-factor variation may be as good 
as 10 %. Note, however, that here we neglected variation 
of the mixing parameter b which is hard to estimate. 

Now we can estimate sensitivity of the ratio of the hy- 
perfine transition frequencies to variation of niq/ Aqcd- 
For ^^^Hg and hydrogen we have 



S[A{Hg)/AiH)] 
[A{Hg)/Am 



2.3- 



Sa 



0.02 



S[mq/AQCD] 
[mq/AQco] 



(17) 



3 



Therefore, the measurement of the ratio of Hg and hy- 
drogen hyperfine frequencies is practically insensitive to 
the variation of masses and strong interaction. The re- 
sult of measurement [T^ may be presented as a limit on 
variation of the parameter a = a[mq/AQC£)]~°-°^: 



|4^|<3.6xlO-"/year 



(18) 



Other ratios of hyperfine frequencies are more sensitive 
to niq/KQCD- For ^^^Cs/^''Rb we have 



5[A{Cs)IA{Rb)] 
[A{Cs)IA{Rh)\ a 



[mq/AQCD] 



Therefore, the result of the measurement may be 
presented as a limit on variation of the parameter X — 



a 



0.49 



mq/AQco] 



0.17. 



The optical clock transition energy E{Hg) (A=282 nm) 
in Hg+ ion can be presented in the following form: 



E{Hg) = const x [—^]Frei{Za) 



(23) 



Note that the atomic unit of energy (first bracket) is can- 
celed out in ratios, therefore, we should not consider its 
variation. Numerical calculation of the relative variation 
of E{Hg) has given 0|: 



5E{Hg) 
E{Hg) 



-3.2 



5a 



(24) 



Variation of the ratio of the Cs hyperfine splitting A{Gs) 
to this optical transition energy is equal to 



1 dX 
'X~dt 



(0.2 ± 7) X 10-^72/ear 



(20) 5[A{Cs)/E{Hg)] 
[A{Cs)/E{Hg)] 



Note that if the relation Q is correct, variation of X 
would be dominated by variation of [mq/ Kqqd]- The 
relation (j^J would give X oc and limit on a variation 
if = (0.03 ±1) X lO-ie/year . 
For i33Cs/H we have 



5[A{Cs)/A(H)] 
[A[Cs)/A{H)] 



5a 5\'m„/Kocn\ 
0.83— 0.2 ; ^' ; (21) 

a {mqlkqcD] 



Therefore, the result of the measurements [I^ may be 
presented as a limit on variation of the parameter Xfi = 



,,0.83 r 



, 1 dXa . 
^Xh dt ' 



< 5.5 X 10 ^^/year 



(22) 



If we assume the relation (O, we would have Xh oc a^, 
\1.^\< 0.7 xlO-^y year. 



a [me/AQco] 



[mq/AQCD] 
(25) 



Here we have taken into account that the proton mass 
Mp oc Aqcd- The factor 6.0 before 5a appeared from 
a^Frei in the Cs hyperfine constant (2-1-0.83) and a- 
dependence of E{Hg) (3.2). Therefore, the work 
gives the limit on variation of the parameter U = 



[melI^QCD\[mq/AQCD\ 



15 



' I year 



(26) 



oc a 



If we assume the relation Q, we would have U 

< ^ lO^^^/year. Note that we presented 
such limits on ^ f as an illustration only since they are 
strongly model-dependent. 

This work is supported by the Australian Research 
Council. 



[1] J. K. Webb , V.V. Flambaum, C.W. Churchill, M.J. 
Drinkwater, and J.D. Barrow, M.J. Drinkwater, and J.D. 
Barrow, Phys. Rev. Lett., 82, 884-887, 1999. J.K. Webb, 
M.T. Murphy, V.V. Flambaum, V.A. Dzuba, J.D. Bar- 
row,C.W. Churchill, J.X. Prochaska, and A.M. Wolfe, 
Phys. Rev. Lett. 87, 091301 -1-4 (2001). M. T. Mur- 
phy, J. K. Webb, V. V. Flambaum, V. A. Dzuba, C. 
W. Churchill, J. X. Prochaska, J. D. Barrow and A. 
M. Wolfe, Mon.Not. R. Astron. Soc. 327, 1208 (2001) ; 
|astro-ph/0012419 M.T. Murphy, J.K. Webb, V.V. Flam- 
baum, C.W. Churchill, and J.X. Prochaska. Mon.Not. R. 
Astron. Soc. 327, 1223 (2001); astro-ph/0012420 M.T. 
Murphy, J.K. Webb, V.V. Flambaum, C.W. ChurchiU, 
J.X. Prochaska, and A.M. Wolfe. Mon.Not. R. Astron. 
Soc. 327,1237 (2001); astro-ph/0012421 

[2] J-P. Uzan. hep-ph/0205340 

[3] P. Langacker, G. Segre and M.J. Strassler, Phys. Lett. 



B 528, 121 (2002); hep-ph/0112233 See also X. Cal- 
met and H. Fritzsch, Eur. Phys. J, C24, 639 (2002); 
hep-ph/0112110 W.J. Marciano, Phys. Rev. Lett. 52, 
489 (1984). 

[4] V.V. Flambaum, E.V. Shuryak, Phys. Rev. D65, 103503 

1-11 (2002); hep-ph 0201303 
[5] K.A. Olive, M. Pospelov, Y.-Z. Qiang, A. Coc, M. 

Casse, E. Vangioni-Flam, Phys. Rev. D66, 045022 (2002); 

hep-ph/0205269 
[6] V.F. Dmitriev, V.V. Flambaum, accepted to Phys. Rev. 

D; astro-ph/0209409 
[7] V.V. Flambaum, E.V. Shuryak, submitted to Phys. Rev. 

D; hep-ph 0212403 
[8] M.T. Murphy, J.K. Webb, V.V. Flambaum, M.J. 

Drinkwater, F. Combes and T. Wiklind. Mon.Not. R. 

Astron. Soc. 327, 1244 (2001); astro-ph/0101519 
[9] L.L. Cowie and A. Songalia, Astrophys. J. 453, 596 



4 



(1995). 

[10] A.I.Shlyakhter, Nature 264 (1976) 340; T.Damour 
and F.J.Dyson, Nucl.Phys.B 480 (1996) 37. Y.Fujii, 
A.Iwamoto, T.Fukahori, T. Ohnuki, M. Nakagawa, H. 
Hidaka, Y. Oura, P. Moller. Nucl.Phys.B 573 (2000) 377. 
[11] H. Oberhummer, R. Pichler, A. Csoto, nucl-th/9810057| 
[12] S.R. Beane and M.J. Savage, hep-ph/0206113 
[13] S. G. Karshenboim. Canadian Journal of Physics 78, 639 
(2000); physics/0008051 S. G. Karshenboim. In Laser 
Physics at the Limits, ed. by H. Figger, D. Meschede 
and C. Zimmermann (Springer- Verlag, Berlin, Heidel- 
berg, 2001) pp. 165-176; physics/0201050 
[14] J.D. Prestage, R.L. Tjoelker, and L. Maleki. Phys. Rev. 

Lett 74, 3511 (1995). 
[15] N.A. Demidov, E.M. Ezhov, B.A. Sakharov, B.A. Ul- 



janov, A. Bauch, and B. Fisher, in Proceedings of the 
6th European Frequency and Time Forum. Noordwijk, 
the Netherlands, 1992 (European Space Agency, Noord- 
wijk, 1992), pp. 409-414. L.A. Breakiron, in Proceedings 
of the 25th Annual Precise Time Interval Applications 
and Planning Meeting, NASA conference publication No. 
3267 [U.S. Naval Observatory Time Service Department 
(TSSl), Washington DC, 1993], pp. 401-412. 

[16] V.A. Dzuba, V.V. Flambaum, J.K. Webb. Phys. Rev. A 
59, 230 (1999). V.A. Dzuba, privite communication. 

[17] H. Marion, et al, physics/0212112j 

[18] S. Bize et al, physics/0212109 

[19] D.B.Leinweber, D.H. Lu, and A.W. Thomas, Phys. Rev. 
D60 , 034014 (1999); hep-lat/9810005, 



